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SUMMARY 

Measurements  made  at  a  Mach  number  of  0.18  and  a  chord-based  Reynolds 
number  of  4.2  »  10^  on  a  constant-chord  model  having  a  NACA  4412  aerofoil  section 
are  described  and  compared  with  the  results  of  flow  field  calculations. 

Both  the  experimental  arrangement  and  the  difficulties  initially  experienced 
in  achieving  an  adequate  approximation  to  two-dimensional  flow  above  the  wing  are 
briefly  outlined. 

The  measurements  include  static  pressure  distributions  on  the  wing  surface 
and  on  the  wind  tunnel  walls  above  and  below  the  mid-span  section  of  the  wing. 

The  main  emphasis  in  the  experiment  was,  however,  on  defining  the  development  of 
the  upper  surface  boundary  layer  through  separation  (at  about  20Z  chord  ahead  of 
the  trai ling-edge)  and  on  into  the  wake,  making  extensive  use  of  laser  anemometry. 

The  flow  field  calculations  are  of  the  semi-inverse  kind  in  which  an  inverse 
momentum-integral  treatment  of  the  shear  flow,  used  to  avoid  difficulties  at 
separation,  is  coupled  to  a  direct  solution  of  the  inviscid  flow  problem.  The 
main  .matures  of  the  method  are  outlined. 
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1  Introduction 

It  1*  now  well  established  chat  flows  past  meny 
types  of  bodies  on  which  there  are  attached  bouny 
dary  layers  can  be  computed  In  a  variety  of  ways1 
and  lntenalve  efforts  are  now  being  aeda  to  extend 
computational  capability  to  cover  flows  containing 
reglona  of  flow  separation.  Such  extanalons  will 
clearly  have  conelderable  practical  Impact  on 
aerodynamic  design.  One  example,  with  which  the 
present  work  Is  concerned,  occurs  In  aerofoil  pro¬ 
file  development  where  the  ability  to  predict  the 
stalling  behaviour  of  given  aerofoils  will  be  most 
valuable. 

A  technique  which  has  proved  successful  in  the 
estimation  of  the  effects  of  viscosity  In  the 
attached  flow  around  aerofoils  Is  to  match  a 
calculation  of  the  outer  lnvlscid  flow  with  a 
calculation  of  the  boundary  layer  development  on 
the  surface  of  the  aerofoil.  In  a  recent  paper 
Williams2  has  drawn  attention  to  certain  obstacles 
to  the  extension  of  this  technique  to  separated 
flows  and  illustrated  ways  in  which  they  any  be 
overcome  in  calculations  for  incompressible 
lnvlscid  outer  flow.  Although  most  of  the  Impedi¬ 
ments  are  related  to  the  nature  of  the  calculation 
procedure  and  have  to  be  tackled  by  adopting 
suitable  strategies  in  formulating  computational 
methods ,  the  adequate  modelling  of  both  the 
separated  parts  of  boundary  layers  on  aerofoils 
and  of  their  wakes  also  presents  difficulties. 

Information  obtained  from  experiment  Is  a  pre¬ 
requisite  to  the  solution  of  these  viscous  flow 


problems.  The  first  four  of  the  following 
Sections  are  devoted  to  an  experiment  designed  to 
provide  such  lnformstlon  for  the  particular  case 
of  an  NACA  4412  aerofoil  tested  at  an  Incidence 
close  to  that  for  maximum  lift.  This  provides  an 
example  of  the  type  of  stall  which  develops  with 
Increasing  incidence  by  the  progressive  enlarge¬ 
ment  of  a  region  of  separated  flaw  on  the  rear 
part  of  the  upper  surface  of  the  aerofoil.  An 
additional  reason  for  choosing  the  NACA  4412  sec¬ 
tion  for  the  experiment  was  to  create  an  oppor¬ 
tunity  for  comparing  results  with  other  work  such 
as  that  of  Wadcock  .  Some  comparlaons,  not 
repeated  here,  with  Wadcock 'a  work  for  both  theory 
and  experiment  are  included  in  Refs  2  and  4. 

Ref  4  also  gives  additional  details  of  the  RAH 
experiment. 

In  the  second  part  of  the  paper  a  theoretical 
method  for  computing  flows  containing  regions  of 
separation  la  described  and  its  predictions  are 
compared  with  the  experimental  results. 

2  Experimental  Details 

The  experiment  was  performed  in  the  13ft  *  9ft 
(4  a  X  2.74  m)  Low  Speed  Wind  Tunnel  of  the  Royal 
Aircraft  Establishment  at  Bedford.  This  facility 
is  a  conventional,  closed-circuit,  wind  tunnel 
with  excellent  flow  quality,  operating  at 
atmospheric  pressure.  The  tunnel  Is  large  enough 
to  accommodate  models  of  sufficient  size  to  permit 
detailed  study  of  their  flow  fields. 

In  this  case  the  model  has  a  constant  chord  of 
1  m,  a  NACA  4412  aerofoil  section  and  spans  the 
full  width  (4  m)  of  the  tunnel.  It  Is  mounted 
with  its  quarter-chord  point,  about  which  it 
pivots  In  the  pitch  plane,  on  the  tunnel  centre 
line.  It  Is  extensively  Instrumented  with  static 
prasaure  orifices.  Boundary  layer  and  wake 
measurements  were  made  at  mid-span  where  the  main 
chordwise  set  of  88  pressure  orifices  Is  located. 

When  first  tested  at  an  Incidence  of  approximately 
14*,  the  model  was  found  to  have  strongly  three- 
dimensional  flow  over  the  rear  part  of  Its  upper 
surface;  elsewhere  the  flow  was  very  well-behaved. 
Surface  oil-flow  studies  revealed  that  the  three- 
dimensional  flow  near  the  surface  of  the  wing  was 
composed  of  the  kind  of  cellular  entitles  which 
are  now  recognised  to  be  fairly  common5  In 
separated  flows.  A  short  leading-edge  separation 
bubble  was  also  detected,  separation  and  reattach- 
mant  being  roughly  at  IX  and  1.7%  chord  respec¬ 
tively.  As  the  overriding  consideration  In  this 
experiment  was  to  obtain  laser  anemometer  meaaure- 
ments  In  two-dimensional,  separated,  mean  flow  It 
was  Imperative  to  take  measures  to  control  the 
flow  without  Impeding  optical  access  to  the  aid- 
span  part  of  the  model. 

Further  oil-flow  studies  Indicated  that  the 
natural  flow  on  tha  model  was  very  close  to  being 
two-dimensional  up  to  an  Incidence  of  U*  where 
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separation  was  at  about  90Z  chord*  It  was  found 
that  the  region  of  acceptable  flow  quality  could 
be  extended  to  13*  by  adding  judiciously-placed, 
streaawlse  fences  running  around  the  leading  edge 
of  the  wing  and  continuing  along  the  front  half  of 
its  upper  surface .  The  aost  laportant  fences  were 
a  pair  symmetrically  placed  on  either  side  of  aid- 
spen  and  1.7  chords  apart.  Because  cellular  stall 
patterns  appeared  in  the  spaces  between  these  fen¬ 
ces  and  the  sides  of  the  tunnel,  an  additional 
fence  was  placed  in  each  of  these  spaces  to  split 
and  weaken  these  outboard  cells.  The  additional 
fences  were  1.5  chords  f roe  aid-span. 

It  was  also  considered  possible  that  there  sight 
be  irregularities  in  the  leading-edge  bubble 
separation  and  reattacfaaant  process  that  could  be 
partly  responsible  for  che  developaent  of  chree- 
dlaenslonallty  in  Che  flow  further  downstream.  To 
counter  this  possibility  an  unusually  wide, 
boundary  layer  tripping  band  was  fitted  to  the 
upper  surface.  This  consisted  of  densely  packed, 
0.28  an  dlaaeter  (0.011  in)  ballotlnl  attached  by 
a  thin  layer  of  adhesive.  The  band  extended  30  tea 
along  che  aerofoil  contour  froa  che  leading-edge; 
le  to  1.4Z  chord  along  the  chord  line.  Transition 
was  fixed  on  the  lower  surface  by  a  similar,  but 
auch  narrower,  band  extending  froa  10X  to  lit  of 
the  chord. 

A  geoaetrlc  incidence  of  12.15*  was  finally 
salacted  for  che  aaln  series  of  testa  so  that  the 
flow  would  be  steadier  than  at  13*,  where  con¬ 
tinuous  visual  observation  of  developing  oil-flow 
patterns  had  indicated  that  the  flow  over  the  rear 
part  of  the  wing  contained  complicated  transient 
phenomena.  The  undisturbed  free  stream  speed  for 
the  tests  was  61  a/s,  streaa  teaperaturea  (which 
are  not  controllable)  ware  22  ±10*  C,  while  the 
Mach  and  Reynolds  numbers  were  0.18  and  4.2  * 

106  respectively.  Apart  froa  the  main  tests, 
there  were  subsidiary  tests  at  lower  incidences  to 
define  the  aerofoil's  lift  curve. 

The  aaln  instrument  used  for  studying  Cha  boundary 
layers  and  wakes  was  a  two-coaponent,  photon- 
correlation,  laser  Doppler  aneaoaeter  (LDA).  This 
was  operated  in  back-scatter  and  was  mounted  on  a 
traverse  aechanlaa  situated  Just  outside  the  tun¬ 
nel  and  to  starboard  of  the  aodel.  Frequency- 
shifting  by  Bragg  calls  was  Incorporated  to 
resolve  directional  aablgultles  in  aaasured  velo¬ 
cities.  The  two  output  signals  for  each  aeesure- 
aent  ware  fed  first  to  separate  autocorrelators 
and  thence  to  a  digital  computer  for  storage  and 
processing.  The  end-products  were  probability 
distributions  for  the  aagnltudes  of  the  two, 
always  orthogonal,  velocity  components  that  ware 
being  aaasured. 

The  aneaoaeter  is  designed  so  that  the  orientation 
of  the  pairs  of  velocity  coaponents  can  be  varied. 
Coaponents  inclined  both  at  0*,  90*  and  ±45*  to 
the  line  of  traverse  were  always  aaasured.  In 
principle  this  should  enable  the  shear  and  noraal 
coaponents  of  Reynolds  stress  in  the  plana  of  the 
aean,  two-dimensional  flow  to  be  evaluated;  in 
practice  difficulties  were  encountered  which  are 
attributable  ultimately  to  the  remoteness  of  the 
aneaoaeter  froa  the  aaasurlng  plana  coupled  with 
various  practical  constraints  on  the  geoaatry  of 
the  optical  spates'*.  Consequently  only  aean  velo¬ 
city  data  froa  the  LDA  are  preeented  here.  The 


measuring  region  yielding  each  data  point  has  an 
effective  dlaaeter  of  2  a  (less  than  1/30  of  the 
thickness  of  the  separated  boundary  layer)  in  the 
plane  of  the  aean  flow  and  extends  about  20  a a 
spanwlse.  Here  again  the  size  of  the  measuring 
region  le  related  to  its  distance  froa  the 
aneaoaeter. 

The  boundary  layer  and  wake  were  always  traversed 
along  lines  norael  to  the  wing  surface  and  to  the 
(extended)  wing  chord  line,  respectively.  The  air 
streaa  was  seeded  with  atoalsed  oil  droplets 
Introduced  at  the  downstream  end  of  the  working 
section  of  the  tunnel. 

To  coapensate  as  far  as  possible  for  the  lack  of 
Reynolds  stress  data  froa  the  laser  anemometer, 
some  additional  traverses  with  hot-wire  probes 
were  aade.  These  probes  were  mounted  or  a  long 
ara  projecting  forward  froa  a  precision  traverse 
gear6  originally  built  some  years  earlier  for 
other  work  on  wakes.  The  same  equipment  was  used 
tc  obtain  velocity  profiles  froa  pitot  tube  sur¬ 
veys  at  201  and  40Z  chord  on  the  upper  surface  and 
at  99. 7Z  chord  on  the  lower  surface. 

3  Wind  Tunnel  Boundary  Conditions 
and  Interference 

Table  1  lists  static  pressures  aeasured  on  the 
longitudinal  centre-lines  of  the  roof  and  floor  of 
the  working  section.  These  pressures  have  been 
used  in  conjunction  with  the  aethod  of  Ashlll  and 
Weeks7  to  find  the  tunnel-wall-lnduced  velocity 
perturbations  midway  between  the  tunnel  roof  and 
floor  for  plana  flow.  As  well  as  an  average 
blockage-induced,  velocity  increment  of  0.00730^ 
(for  an  undisturbed  flow  speed  of  D^)  there  exist 
both  a  aean  Induced  inclination  of  the  streaa  of 
about  0.34*  and  a  curvature  equivalent  to  a  camber 
of  -0.28Z.  Fig  1  shows  the  induced  upwash 
and  streamline  ordinates  (z^/c). 

The  boundary  layers  on  the  roof  and  floor  of  the 
eapty  working  section  3.29  a  ahead  of  the  position 
of  the  aodel 's  leading  edge  have  the  following 
characteristics:  displacement  thickness  (5*)  - 
7.6  ib,  Reynolds  number  for  momentum  thickness 
0  *  24400,  shape  parazMter  H  (H  •  6*/0)  -  1.30. 

4  Measured  Wing  Static 
Pressure  Distribution 

The  pressure  distribution  around  the  aerofoil  pro¬ 
file  at  aid-span  is  shown  in  Fig  18  where  it  is 
coapared  with  the  results  of  calculations  to  be 
described  in  Sections  6  and  7.  The  measurements 
do  not  define  the  pressure  distribution  just 
behind  the  leading  edge  well  enough  to  yield 
extreatly  accurate  force  coefficients  at  high 
Incidences.  It  is  believed  that  the  lift  coef¬ 
ficient  le  Ci  •  1.46±0.01  at  the  Incidence  of 
12.15*  for  the  aeln  series  of  tests. 

At  this  incidence,  pressures  were  very  accurately 
repeatable  over  the  whole  of  the  lower  surface  of 
the  wing,  where  the  pressure  gradients  were  all 
favourable,  and  also  if  x/c  was  less  than  0.45  on 
the  upper  surface  (x  being  aeasured  along  the 
chord  line  of  length  c).  Further  downstreaa  on 
Cha  upper  surface  differences  between  repeated 
aeasuraaents  begin  to  appear.  These  differences 
asy  be  attributed  to  aoae  unsteadiness  in  the 
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flow,  for  which  Cher*  was  visual  evidence  as  noced 
In  Section  2.  Just  ahead  of  separation  the 
scatter  indicated  by  four  repeated  measurements  la 
0.025  in  C..  At  the  saae  station,  x/c  *  0.75,  the 
static  pressure  data  froa  a  spanwlse  row  of 
measuring  orifices  fell  within  a  band  whose  width 
Is  0.06  in  C_  for  stations  within  ±0.6  c  of  aid- 
span.  This  figure  Includes  both  scatter  and  saall 
systeaatlc  spanwlse  variations.  Although  the 
chordwlse  pressure  distribution  is  not  tabulated 
In  detail,  values  of  C_  at  boundary  layer  traverse 
stations  are  Included  In  Table  2. 

5  Boundary  Layer  and  Wake  Measurements 

Fig  2  gives  a  general  lapression  of  the  scope  of 
the  boundary  layer  and  wake  velocity  profile 
aeasureaents.  Leading  characteristics  of  the  pro¬ 
files  are  listed  In  Tables  2  and  3  for  the  boun¬ 
dary  layers  and  wake  In  turn.  As  aay  be  seen  froa 
Fig  2,  boundary  layer  separation  was  found  to  be 
at  approximately  x/c  ■  0.8  on  the  upper  surface. 
Fig  3  shows  four  aesn  velocity  profiles  obtained 
by  laser  aneaoaetry  In  the  separated  boundary 
layer  on  two  different  occasions  and  with  two 
laser  beaa  orientations.  The  later  aeasureaents 
are  racher  more  consistent  because  the  aneaoueter 
has  been  progressively  Improved  during  the  course 
of  the  experlaent.  Normalised  profiles  for 
0.66  <  x/c  <  0.997  are  shown  In  Fig  A.  The  velo¬ 
city  component  parallel  to  the  wing  surface  (0)  is 
given  as  a  fraction  of  Ue  its  value  at  the  edge 
of  the  boundary  layer;  distance  (n)  normal  to  the 
wing  surface  Is  given  as  a  fraction  of  5,  Its 
value  when  0  •  0.995  0  ;  in  the  Interests  of 
clarity  smoothed  profiles  and  not  Individual  data 
points  are  shown.  The  shapes  of  the  profiles 
change  In  an  orderly  way  as  separation  Is 
approached  and  passed.  It  Is  worth  noting  that 
there  are  measurable  velocity  gradients  normal  to 
the  wing  surface  In  the  lnvlsdd  flow  outside  the 
boundary  layer. 

The  presence  of  these  velocity  gradients  Is  an 
Indication  that  Che  boundary  layers  are  too  chick 
for  classical  first-order  boundary  layer  theory  to 
apply  and  highlights  Che  Importance  of  taking 
account  of  higher  order  effects  In  the  calcula¬ 
tions  which  are  described  In  Sections  6  snd  7. 
There  la  the  further  Implication  that  generallsad 
definitions  of  bound|ry  layer  displacement  and 
momentum  thicknesses9  should  be  used  In  the  pre¬ 
sentation  of  experimental  data.  It  should  be 
noted  that  this  has  not  yet  been  done  for  the  data 
of  this  experlaent;  le  the  quantities  o.  and  0  of 
the  equivalent  lnvlsdd  flow  (EIF)  of  S<ctlons16.1 
and  6.2  have  been  taken  to  be  respectively  a 
coostanc  (for  this  low  Mach  number  flow)  and  0  In 
analysing  the  data.  * 

It  aay  be  seen  froa  Fig  3  that  the  mean  velocity 
profiles  are  subject  to  some  uncertainty,  par¬ 
ticularly  near  the  wing  surface;  this  affects  the 
derived  valuee  Included  in  the  Tables.  Momentum 
thickness  Is  the  quantity  most  seriously  affected. 
Examination  of  the  data  suggests  that  S,  &*,  8  and 
B  are  subject  to  the  following  maximum  uncertain¬ 
ties:  ±5Z,  ±3X,  ±8X  and  ±8X  respectively. 

The  aass  flow  shape  parameter  fl,  ■  («-6»)/8]  can 
also  be  determined.  Values  of  H:  from  this 
experiment  ere  In  satisfactory  agreement  with  the 
curve  of  Fig  14  for  attached  flow;  for  separated 


flow,  however,  smoothed  values  of  H,  appear  to 
remain  constant  at  H,  »  A  when  H  exceeds  5.  In 
order  to  reproduce  the  rising  trend  In  Fig  1A, 
values  of  <  that  gradually  exceeded  the 
experimentally-derived  ones  would  have  to  be  cho¬ 
sen;  for  example  the  required  Increment  would  be 
10X  at  H  *  7.5.  With  regard  to  this  observation. 

It  la  worth  noting  two  points.  Firstly  the  curve 
passes  through  H  ■  A  at  H  ■  A  which  Is  consistent 
with  a  cosinusoidal  wake-llke  profile  If  6  Is 
defined  to  be  where  0  •  U  .  This  version  of  6 
exceeds  6  for  U  -  0.995U  eby  almost  5Z  for  that 
profile.  Secondly  we  have  noticed  a  consistent 
tendency  for  U,  If  It  Is  very  close  to  0  ,  to 
approach  0e  rather  more  quickly  when,  asehere,  It 
Is  measured  by  a  laser  anemometer  than  when  It  Is 
measured  by  an  (Intrusive)  hot-wire  probe  and  so 
the  laser-derived  boundary  layer  thicknesses  may  be 
slightly  less  Chan  those  that  other  instruments 
would  yield. 

Skin  friction  coefficients  for  the  attached-flow 
boundary  layers  have  been  estimated  by  the  Clauser 
chart  method  (for  a  Karman  constant  of  0.A1  and  an 
additive  constant  of  5.2  In  the  logarithmic  part 
of  the  velocity  profile),  though  only  within 
rather  broad  limits  for  flows  near  separation 
where  the  logarithmic  parts  of  the  profiles  are 
very  small.  If  we  denote  the  skin  friction  coef¬ 
ficient  and  the  value  predicted  by  the  well-known 
Ludwleg  and  Tlllmann  formula  by  Cf  and  c^ 
respectively,  the  results  are  well  represented 
by  the  Interpolation  formula  <•  •  0.985  CjL  - 

3  X  10-5.  This  relationship  Is  consistent  with  a 
value  of  A  for  H  at  separation  In  this  experiment. 
The  estimates  of  cf  are  included  In  Table  2. 

Reynolds  stressss  In  the  outer  parts  of  the  boun¬ 
dary  layars  were  measured  with  a  hot-wire  anemo¬ 
meter.  Comparison  of  mean  velocity  profiles 
obtained  by  the  wire  and  laser  anesiometers  showed 
that  there  were  probe  Interference  effects  on  the 
former  and  so  the  Reynolds  stresses  are  presented 
In  Figs  5  and  6  only  for  Che  parts  of  the  boundary 
layers  where  the  two  types  of  anemometer  yielded 
mean  velocity  profiles  of  very  similar  shapes. 

""7  2  2  2 

The  dimensionless  normal  stresses  u  / U  and  w  /D^ 
progressively  increase  In  Importance  from  x/c  " 

0.59  to  the  tralling-edge  of  the  wing.  In 
contrast  the  peaks  of  the  dimensionless  shear 
stressss  Initially  Increase  but  soon  reach  a  maxi¬ 
mum,  at  which  they  remain  from  at  least  x/c  *  0.66 
to  x/c  ■  0.85,  before  falling  slightly  near  the 
trailing  edge.  Peak  shear  stress  coefficients, 

(cT)>>x  »  f-2uw/0x)max,  sre  Included  in  Table  2. 

Since  the  normal  stresses  are  clearly  of 
Increasing  Importance  to  the  overall  momentum 
balance  of  the  boundary  layer  near  and  beyond 
separation,  an  attempt,  recorded  In  the  final 
column  of  Table  2,  has  been  made  to  estimate  their 
contribution  to  the  momentum  Integral  equation. 

The  estimate  relies  on  the  sssumptlons  that  the 

she  pea  of  the  distributions  of  -  w^  )/D*  sre 

mainly  functions  of  (H-R  )/H  and  may  be  scaled  by 
o 

using  the  maximum  values  of  fu^  -  H  is 

the  value  of  H  In  a  datum  equilibrium  flow  with 
aero  pressure  gradient  at  the  relevant  local  value 
of  Rg.  The  tabulated  quantity  1/8  is  defined  by 
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ics  algebraic  significance  la  that  one  nay  regard 
it  as  reducing  8  to  8(1  •  1/6)  in  the  conventional 
*o men tun  Integral  equation.  This  nay  be  written1, 
with  normal  stresses  (but  no  other  higher-order 
terms)  Included,  as 

h  fo^fe-l)]  +  «*3;(*0*)  -  ictf  ; 


here  I  la  the  nonaal  atress  Integral  defined 
above.  Integrating  fro.  an  Initial  atatlon  Xj 

x,  and  dividing  by  the  value  of  9  at  x, ,  we 
obtain  e 


0^  9-1) 


2f0f9)l 


i: 


cfof-J 


dU 

_ t 

dx 


dx 


If  I  la  negligible  at  x^. 


to 


A  momentum  balance  for  the  upper  aurface  boundary 
layer  la  given  In  Fig  7.  A  datum  at  x/c  >0.2  has 
been  choaen.  The  poaltlve  contribution  la 

uf8/(uf9)j-l,  ahovn  by  aolld  line*.  The  effect  of 


■edifying  9  to  take  account  of  normal  stresses  aa 
juat  prevloualy  outlined  la  Indicated  by  broken 
lines.  The  balancing  contrlbutione  (below  the  x/c 
axis)  coepri.se  l~ 

d  erf  | 


1 


!x 

0.2c 


- 


4*  . 


dx 


dx. 


The  flow  la  seen  to  be  nearly  In  balance  In 
general  and  so  ahould  represent  a  good  approxlma- 
tlon  to  two-dimensional  mean  flow;  however  the 
scatter  of  the  values  of  9  and  4*  makes  It  dif¬ 
ficult  to  draw  firm  conclusions-  Beyond  separa¬ 
tion  (x/c  •  0.8)  the  contribution  of  the  estimated 
normal  stresses  la  significant  and  It  would 
clearly  be  worthwhile  to  Incorporate  estimates  of 
other  higher  order  effects. 


because  the  velocity  gradients  are  very  steep  and 
so  a  very  slight  change  In  the  position  of  the 
streamline  has  a  disproportionately  large  effect 
on  the  estimated  velocity.  Further  downstream  the 
velocity  on  the  dividing  streamline  Increases  as 
mixing  begins  to  equalise  the  velocities  within 
the  wake. 

Since  the  upper  and  lower  part*  of  the  wake  are 
derived  respectively  from  a  very  thick  separated 
boundary  layer  and  a  very  thin  attached  boundary 
layer  they  are  Initially  vary  different  and  both 
change  markedly  just  behind  the  wing.  Reversed 
mean  flow  persists  In  the  upper  part  of  the  wake 
until  x/c  reaches  1.09.  Aa  Fig  11  shows,  the 
displacement  thlcknese  (4*)  of  the  upper  wake 
falls  and  that  (4*)  of  thX  lower  wake  rises  In 
this  region.  The1 corresponding  momentum 
thicknesses  8  and  9g  do  not  change  as  smoothly. 
The  former  haX  a  maximum  near  x/c  ■  1.2.  This  la 
tsntatlvely  associated  with  the  closure  of  the 
bubble  of  reverse  flow.  This  Interpretation  la 
consistent  with  results  obtained  by  Green’  in  an 
attempt  to  modal  a  flow  that  Is  qualitatively 
similar  In  some  respects  -  namely  reattachnent 
behind  a  rearward-facing  step.  He  found  that  8 
peaked  In  his  calculations  at  a  point  of  maximum 
static  pressure  Just  beyond  reattachment.  The 
peak  found  in  this  experiment  Is  likewise  beyond 
bubble  eloaure  and  roughly  where  the  near-wake 
pressure  recovery  ends  (Figs  2  and  10). 

The  displacement  surfaces  above  the  rear  part  of 
the  wing  and  in  the  wake  are  Included  In  Fig  2. 
Despite  the  rapid  variations  in  displacement 
thickness  near  the  cralllng  edge,  the  upper 
displacement  surface  Is  smooth  and  only  slightly 
curved.  The  net  affect  of  the  displacement  sur¬ 
faces  on  the  outer  lnvlsdd  flow  la  equivalent  to 
a  reduction  In  camber  and  In  Increase  In  thickness 
near  the  trailing  edge.  The  wake  Is  noticeably 
Inclined  to  the  wing  chord  line  within  the  range 
of  the  measurements.  Table  3  Includes  the  esti¬ 
mated  position  of  the  dividing  streamline  in  the 
wake. 


The  development  of  the  upper  surface  boundary 
layer  Is  Illustrated  In  quantitative  terms  In 
Figs  16,  20  and  21.  The  displacement  thickness 
grows  steadily  along  the  chord  and  does  so  par¬ 
ticularly  rapidly  near  and  beyond  separation 
(Fig  20).  This  leads  to  a  rapid  Increase  In  H 
(Fig  21)  over  the  rear  part  of  the  chord  because  9 
(Fig  16)  varies  comparatively  slowly  aft  of 
separation. 

The  wake  has  been  divided  Into  upper  and  lower 
parts  by  estimating  the  position  of  the  mean 
streamline  separating  the  flows  that  have  passed 
above  and  below  the  wing.  This  was  done  by  con¬ 
sidering  the  mess  flux  Into  and  within  the  lower 
and  thinner  part  of  the  wake  on  the  assumption 
that  the  mean  flow  was  two-dimensional. 
Dimensionless  velocity  profiles  for  both  parts  of 
the  wake,  similar  to  those  given  In  Tig  6  for  the 
wing  boundary  layer,  are  shown  In  Figs  8  and  9. 
Dimensionless  velocities  at  the  edges  of  the  wake 
and  on  Its  Internal  dividing  streamline  are  given 
In  Fig  10. 

It  can  be  seen  from  Fig  10  that  there  Is  some 
pressure  recovery  In  the  near  wake,  to  slightly 
beyond  the  value  In  the  undisturbed  stream.  The 
velocities  on  the  dividing  streamline  at  stations 
close  to  the  wing  are  scattered  and  unreliable 


6  General  Description  of  the  Prediction 
Method  for  Separated  Flow 

A  vlseous-lnvlscld  Interaction  method  Is  used.  It 
differs  from  the  type  that  la  successful  In 
treating  attached  flow,  for  which  It  ha*  been 
shown  that  the  effects  of  viscosity  can  be  esti¬ 
mated  by  coupling  a  solution  of  the  outer  lnvlscld 
flow  with  a  direct  solution  for  the  boundary 
layers  and  wake*  associated  with  an  aerofoil.  For 
a  direct  solution  of  the  boundary  layer,  the 
pressure  distribution  Is  specified  and  the  displa¬ 
cement  thickness  (among  other  quantities)  Is 
calculated.  As  the  boundary  layer  approaches 
separation  the  equations  describing  It  become 
stiff10  In  the  sense  chat  a  small  change  In 
pressure  distribution  leads  to  a  large  change  In 
displacement  thickness.  This  stiffness  leads  to 
numerical  problems  In  Che  integration  of  the 
equations,  thus  limiting  the  approach  described  In 
kaf  1  to  attached  flows.  If  the  boundary-layer 
equation*  are  solved  In  invars*  form  -  le  the 
pressure  distribution  Is  calculated  from,  say,  the 
displacement  thickness,  then  close  Co  separation 
the  equations  are  no  longer  stiff  and  Integration 
of  the  equations  presents  no  numerical  dif¬ 
ficulties.  An  Integral  formulation  of  an  Inverse 
turbulent-boundary  layer  method  Is  given  by  East 
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et  al 1 1  and  la  derived  fro*  Che  direct  lag- 
entralnaent  see hod  given  by  Green  et  al12:  an 
extenalon  of  thia  Inverte  aethod  la  uaed  In  the 
present  calculation. 

The  Inverse  solution  of  the  boundary  layer 
equations  Is  aatched  to  a  direct  calculation  of 
Che  lnvlscld  flow  by  a  seal-inverse  scheae  of  the 
sort  described  by  Carter13  and  be  Balleur14.  A 
full  description  of  the  aethod  Is  given  In  Refs  2 
and  IS.  but  an  outline  of  the  aethod  Is  given  In 
the  following  sections.  First  a  description  of 
the  boundary  conditions  and  the  calculations  for 
the  outer  lnvlscld  flow  Is  given.  This  Is 
followed  by  a  description  of  the  boundary  layer 
equations  including  soae  second  order  terns  and 
the  aodif lcatlons  needed  for  separated  flow. 
Finally  a  description  of  the  seal-inverse  scheae 
and  a  relaxation  foraula  derived  froa  a  stability 
analysis  art.'  presented. 


where  <C*  Is  the  curvature  of  the  dlsplaceaent  sur¬ 
face  and  6  Is  the  aoaentua  thickness  of  the  boun¬ 
dary  layer.  Thus  the  pressure  calculated  in  the 
EIF  has  to  be  corrected  before  coaparlson  with 
experlaental  values. 


In  the  wake  slallar  boundary  conditions  apply  on  a 
convenient  line  dividing  the  flows  from  the  upper 
and  lower  surfaces  of  the  aerofoil.  It  Is  helpful 
In  soae  ways  to  regard  this  line  as  a  vortex 
sheet  in  the  EIF  and  refer  to  It  as  such.  The 
dlsplaceaent  effect  of  the  wake  Is  given  by  a  Jump 
aw.  In  the  component  of  velocity  normal  to  the 
wake,  which  by  analogy  with  equation  (1)  Is  given 
by 


AW, 


lw 


U  «*) 

iw  t' 


(3) 


where  suffix  t  denotes  the  total  for  the  wake 
and  ”  denotes  aean  values. 


6.1  Boundary  Conditions  for  the  Outer  lnvlscld 
Flow 

Ref  1  describes  the  matching  of  flows  Including 
che  effects  of  higher-order  terms;  these  terms 
become  more  Important  in  separating  and  separated 
flow.  In  the  matched  solution  It  Is  assumed  chat 
the  flow  field  can  be  split  Into  two  regions;  an 
Inner  region  described  by  the  Navler-Stokes 
equations  (or  some  approximation  of  these 
equations)  and  an  outer  region  described  by  the 
Euler  equations,  which  for  lrrotatlonal 
Incompressible  flow  are  approximated  by  Laplace's 
equation.  If  an  Integral  formulation  Is  used  to 
describe  Che  development  of  Che  boundary  layer 
then  It  Is  most  convenient  to  match  the  inner  and 
outer  solutions  on  Che  aerofoil  surface  and  some 
convenient  line  In  the  wake:  the  definition  of  the 
lnvlscld  flow  Is  extended  through  the  boundary 
layer  and  wake  to  this  matching  line.  Following 
Lock  and  Flrmin  the  new  lnvlscld  flow  la  called 
che  equivalent  lnvlscld  flow  (EIF)  whilst  the  name 
real  viscous  flow  (RVF)  Is  reserved  for  the  physi¬ 
cal  flow.  The  boundary  conditions  for  the  EIF  are 
obtained  by  differencing  continuity  and  momentum 
equations  for  the  EIF  and  RVF,  (the  'deficit 
formulation')  and  are  summarised  In  Fig  12.  The 
transpiration  velocity  Wll#  Is  given  exactly  by 
differencing  the  continuity  equation  for  the  RVF 
and  EIF,  then  Integrating  the  difference  across 
che  shear  layer,  to  give 

w,  •  — 4-  (p,  u,  «*)  (1) 

lw  ds  1  iw  lw  ‘ 

with8  4*  •  L  /"  (o.O,  -  oU)dn  and  s,  n  as 

°lw  lw 

distance  along  and  perpendicular  to  the  surface, 
where  suffices  1  and  w  Indicate  the  equivalent 
lnvlscld  flow  and  conditions  at  the  surface  of  the 
aerofoil  respectively.  Other  symbols  take  the 
usual  aaaning,  0,  W  as  velocity  along  s  and  n,  p 
as  density  and  4*  as  a  displacement  thickness. 


By  neglecting  small  terms  In  the  Integrated  form 
of  the  difference  of  Che  normal  aoaentua  equation 
for  the  EIF  and  RVF,  a  simple  numerical  approxima¬ 
tion  for  the  difference  in  wall  pressure  between 
the  EIF  and  RVF  la  derived  as9: 


Plw-  P» 


*1-  °lw 


K*(«  ♦  «») 


(2) 


The  pressure  variations  across  the  wake  in  the  RVF 
and  EIF  are  different.  The  flows  outside  the  wake 
are  aatched  by  placing  a  vortex  sheet  In  the  wake 
In  the  case  of  the  EIF.  By  Integrating  the  dif¬ 
ference  forms  of  the  normal  momentum  equation  for 
the  EIF  and  RVF  across  the  wake,  the  pressure  dif¬ 
ference  across  the  vortex  sheet  Is  found  to  be: 


where  •  denotes  mean  values  and  Is  the  cur¬ 
vature  of  the  dividing  streamline. 

6.2  Outer  Calculation:  Equivalent  lnvlscld  Flow 

Although  the  matching  scheme  will  allow  solutions 
for  compressible  flow,  for  the  present  low-speed 
calculations  about  a  NACA  4412  aerofoil  the  outer 
flow  Is  assumed  to  be  Incompressible  and  irrota- 
tlonal  and  la  given  by  a  solution  of  Laplace's 
equation.  The  outer-flow  problem  is  now  linear 
and  a  solution  of  the  EIF  Is  obtained  without 
Iteration.  Through  Green's  third  Identity  the 
solution  of  Laplace's  equation  can  be  expressed  as 
a  combination  of  fundamental  solutions  distributed 
on  the  surface  of  the  aerofoil.  In  two  dimensions 
the  fundamental  solutions  are  taken  as  sources  and 
vortices  and  are  approximated  by  piecewise  con¬ 
tinuous  functions  defined  on  the  N  sides  of  a 
polygon  Inscribed  In  the  aerofoil  (see  Fig  13)., (A 
numerically  stable  scheme,  developed  by  Newling1  , 
has  the  source  distribution  as  constant  on  each 
panel  and  the  vortlclty  varying  linearly.  The 
particular  mix  of  sources  and  vortices  Is  chosen 
so  that  the  singularity  densities  on  opposite 
panels  on  Che  upper  and  lower  surface  are  equal. 
The  number  of  unknowns  is  N+l  of  which  N  are 
determined  by  satisfying  a  Neumann  boundary  con¬ 
dition  at  che  centre  of  each  panel,  that  the  nor¬ 
mal  velodcy  should  be  given  by  equation  (1).  The 
final  unknown  Is  Che  strength  of  the  vortlclty  at 
the  trailing  edge  and  smooch  flow  at  the  trailing 
edge  Is  obtained  by  setting  the  vortlclty  on  che 
aerofoil  at  Che  trailing  edge  equal  to  the  vor¬ 
tlclty  obtained  by  approaching  the  trailing  edge 
from  upstream  through  the  wake:  the  value  Is 
obtained  froa  equation  (4). 

The  shape  of  the  streamline  In  che  EIF  froa  the 
trailing  edge  la  calculated  by  an  Interactive  pro¬ 
cess.  An  initial  guess  at  the  shape  of  the 
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streamline  1*  split:  Into  panels  and,  by  sweeping 
down  the  wake,  each  panel  is  progressively  aligned 
with  the  local  flov  direction.  In  the  wake  the 
boundary  conditions  for  the  EIF  are  represented  by 
constant  source  panels  along  the  dividing 
streamline  for  the  displacement  effect  and 
linearly  varying  vortlcity  for  the  curvature 
effects.  The  strengths  of  the  sources  and  vor¬ 
tices  are  derived  from  equations  (3)  and  (4) 
respectively.  In  the  present  method  It  has  not 
been  possible  to  obtain  converged  solutions  with  a 
representation  of  tha  curvature  effects  In  the 
wake,  so  all  Che  calculations  reportsd  have  the 
vortlcity  In  the  wake  set  to  taro. 


6.3  Inner  Calculation:  Shear  layer 


A  set  of  Integral  equations  for  calculating  the 
development  of  a  turbulent  boundary  layer  Is  given 
by  Lock1.  These  equations  use  a  normalised 
transpiration  velocity  S  »  V  PilfD.wd/ds(PiwU^HJ*> 
as  Independent  variable  as  opposed  to  the  displa¬ 
cement  thickness  used  by  East  et  al11.  Both 
Integral  formulations  require  a  correlation 
between  the  shape  factors 


H 


1 
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p-)dn  and  Hj 
lw 


In  Che  Integrals  above,  Che  upper  limit  6  Is  the 
true  edge  of  the  boundary  layer,  not  where 
D  -  0.9950  as  In  the  presentation  of  the  experi¬ 
mental  results  In  earlier  Sections. 


The  new  sec  of  equations  based  on  S  permits  the 
use  of  a  correlation  which  has  a  minimum  close  to 
separation  and  which  consequently  accords  more 
closely  with  experimental  measurements  of 
separating  flow.  The  new  equations  als^lnclude 
higher  order  terns  as  proposed  by  Weeks  The 
four  equations  are  given  below. 

Streanelae  Momentum  Integral  Equation 
.D  C,  ,  a  dU,  ,  ,  9 


d0  .  f  /HW-O-M2: 
ds  1  (H+2'M  > 


— „2  5-{p,  0?  f 0} 
plwUiwd»  lw  lw 


where  c.  Is  the  skin-friction  coefficient,  M  Is 
the  local  Mach  number;  f  embodies  the  second  order 
terms  as  f  ■  K*(9+«*)/28  -  (I/9)ct/ctKq  with  1  »* 

the  Reynolds  normal  stress  Integral  /d[(u2-  w2]/ 

U?ldn,  C  the  shear  stress  coefficient  and  CT_. 
Its  value'ln  equilibrium  conditions.  The  equation 
can  be  reduced  to  standard  form  by  defining 


Entrainment  Equation 

D9H-  V  -tm2V?  - 

where  Cg  Is  the  entrainment  coefficient  and 
D  -  Hj  -  H  dHj/dff  +  (Hj+2dHj/dH 1m2 /5. 

Ml  -entrainment  Equation 


“3.  rfCE,cf 
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where  c.  1*  the  ekln  friction  coefficient  In 
equlllbriua  flov  In  eero  pressure  gradient  and  X 
is  a  scaling  factor  on  the  dissipation  length. 
The  suffix  EQ  denotes  equilibrium  conditions  and 
EQq  equilibrium  In  the  absence  of  secondary 
Influences  on  turbulence  structure. 


Velocity 


0  dDiw  dHl  I  17  dH, 

(H+1)d| - -in  -  — i  S  +  i  c*  fl4d2]H.-H  — i 

Uiw  d*  dH  2  f  L  5  1  dH  J 
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The  extraneous  Influences  on  the  turbulence  struc¬ 
ture  ere  contained  In  the  scaling  factor  X  In  the 
lag  equation  and  the  most  Important  influence  in 
this  context  Is  the  effect  of  longitudinal  cur¬ 


vature  on  turbulence  structure  for  which  X  Is  set 
to  1  +  B^R, (l+M  )/5.  R1  is  the  'Richardson  num¬ 
ber’  for  which  an  average  value  through  an  equili¬ 
brium  boundary  layer  la  given  approximately  by 


R1  -  ilfH+HjK  -!  +  0.3]  ; 


B.  la  taken  from  Ref  12  as  7  on  a  convex  wall  and 
475  on  a  concave  wall  while  R  la  the  radius  of 
longitudinal  curvature. 

The  method  Is  completed  by  the  specif lcatlon  of  a 
correlation  between  H  and  H  which  can  be  derived 
from  velocity  profiles.  For  the  calculation  of 
separated  flow  velocity  profiles  must  contain 
some  reverse  flow.  The  two  parameter  jo»  speed 
velocity-profile  family  of  Le  Balleur  and  Lock1 
provides  a  suitable  profile  and  la  represented 
by  the  formula 

§•  -  1  +  Cjlogn  -  C2F(n) 

e 

where  n  »  y/S  and  F(n)  takes  the  form 
F(n)  »  0  If  0  <  n  <  r\*  , 

F(n)  •  ff-f~^)  If  n*  <  n  <  1  , 

and  n*  is  rero  for  attached  flows  but  defined 
empirically  for  separated  flows.  Lock  found  chat 
a  close  fit  to  experimental  profiles  is  obtained 
if  f  Is  taken  as  Coles'  wake  function,  so  that 

f(0  ■  Kl+cossC) 
where  E  ■  (n-n*)/(I-n*). 

These  profiles  can  be  used  to  determine  a  rela¬ 
tionship  betwejg  H  and  Kj  which  Is  Illustrated  In 
Fig  14.  Green  _hes  argued  that  the  'transformed' 
shape  parameter  H  might  be  regarded  as  the  equiva¬ 
lent  In  compressible  flows  of  the  parameter  H  at 
low  speeds,  thus  this  (H,  H, )  relationship  can  be 
used  for  compressible  flows; 


6.4  Vlscous-lnvlscld  Matching 


A  flow  diagram  for  tha  semi-inverse  method  of 
matching  the  inner  and  outer  flows  Is  given  In  Fig 
15.  An  Initial  estimate  of  the  normalised 
transpiration  velocity  S  is  used  In  an  inverse 
calculation  of  the  boundary  layer  to  give  an  esti¬ 
mate  of  the  velocity  distribution  In  the  EIF. 
The  same  estimate  Is  also  used  as  *  boundary  con¬ 
dition  in  the  direct  lnvlscld  method  to  yield 
another  estimate  U1  of  the  velocity  In  tha  EIF. 
The  difference  between  the  two  estimates  of  tbs 


V  I 

velocity,  -  U  ,  is  used  to  correct  the  nor¬ 
malised  transpiration  velocity  and  the  two 
simultaneous  calculations  are  repeated*  A  stabi¬ 
lity  anal^gi^of  a  linearised  form  of  the 
equations'  »  which  examines  the  growth  of 
errors  in  the  normalised  transpiration  velocity 
and  streamwise  velocity  expressed  as  Fourier  modes 
yields  a  correction  formula  for  S, 


AS  * 


0B8 

"l  dUIw 

i 

duH 

iw 

B+B0./AX 

ds 

"u? 

lw 

ds 

which  proves  useful  in  the  full  non-linear 
problem.  The  factor  6B9/( 8+B0x/Ax)  can  be  viewed 
as  a  relaxation  factor  which  varies  throughout 
the  fleld.of  calculation  and  depends  upon 
8  *  /(1-M  j,  Ax  the  panel  length  in  the  invlscld 
calculation  and  B  -  (H+1)D  dH/dH  ,  ie  quantities 
from  the  boundary  layer. 


6.5  Outline  of  the  Application  of  the  Method 

The  method  must  be  given  the  aerofoil  ordinates, 
the  Reynolds  number,  Mach  number  and  angle  of 
incidence.  Then  the  laminar  portion  of  the  boun¬ 
dary  layer  is  calculated  by  a^ompressible  version 
of  the  method  due  to  Thwaltes  •  The  computer 
program  allows  for  three  different  types  of  tran¬ 
sition  to  a  turbulent  boundary  layer.  Natural 
transition  is  predicted  by  Granville’s 
correlation2  ,  If  the  laminar  boundary  layer 
separates  before  natural  transition  has  been  pre¬ 
dicted  then  the  development  of  the  laminar  separa¬ 
tion  bubble  is  calculj£ed  by  Horton's 
semi -empirical  method2  .  Transition  can  also  be 
fixed  at  a  specified  point  as  long  as  neither  of 
the  other  criteria  has  been  satisfied  upstream. 


Downstream  of  transition  the  development  of  the 
boundary  layer  and  wake  are  calculated  by  the 
Inverse  method  described  in  Section  6.3. 


7  Comparison  Between  Theory  and  Experiment 

Calculations  of  the  flow  about  a  NACA  4412  aero¬ 
foil  are  compared  with  the  experimental  results  at 
a  Reynolds  number  of  4.17  *  10*  and  a  Mach  number 
of  0.18.  In  the  theoretical  calculation  the  tran¬ 
sition  from  laminar  to  turbulent  boundary  layer  is 
fixed  at  x/c  *  0.014  and  0.110  on  the  upper  and 
lover  surfaces  respectively,  which  correspond  to 
the  downstream  ends  of  the  transition  trips  used 
in  the  experiment.  The  theoretical  method  does 
not  implicitly  contain  an  estimate  of  the  change 
in  the  state  of  the  boundary  layer  as  it  passes 
over  the  trip,  but  this  can  be  simulated  by 
increasing  the  momentum  thickness  at  the  trip  on 
the  upper  surface:  at  an  angle  of  incidence  of 
12.15*,  which  corresponds  approximately  to  the 
maximum  lift  coefficient  for  the  aerofoil  in  the 
experimental  tests,  an  increase  in  the  momentum 
thickness  of  the  boundary  layer  of  0.0002  m  at 
x/c  ■  0.014  on  the  upper  surface  leads  to  reaso¬ 
nable  agreement  with  the  momentum  thickness 
measured  at  x/c  *0.2  and  0.4  as  indicated  in  Fig 
16.  On  the  lower  surface  there  is  no  Increase  In 
momentum  thickness  at  the  trip  but  the  calculated 
value  of  0.0008  a  compares  quite  well  with  the 
measured  value  of  0.00064  m,  at  the  trailing 
edge. 


curvature  has  been  Included  but  the  second  order 
effects  of  flow  curvature  and  Reynolds  normal 
stress  have  not.  Calculations  with  the  same 
increase  in  momentum  thickness  at  the  upper  sur¬ 
face  trip  are  repeated  for  angles  of  Incidence 
between  0*  and  14.5*:  the  lift  coefficients  are 
compared  with  experimental  and  lnvlscid  values  In 
Fig  17.  Up  to  an  incidence  of  12.15**  there  is 
good  agreement  with  the  experimental  results. 
However  the  theoretical  results  extend  beyond  the 
range  of  experimental  values  and  indicate  a 
decrease  in  the  lift  coefficient  as  the  separation 
point  moves  towards  the  leading  edge. 

The  calculated  and  measured  pressure  distributions 
at  an  angle  of  incidence  of  12.15*  are  compared  In 
Fig  18.  The  interaction  calculation  was  run  for 
200  Iterations  so  that  it  is  fully  converged  in 
terms  of  lift  coefficient.  The  two  pressure 
distributions  are  similar  except  that  the  experi¬ 
ment  indicates  slightly  lower  pressures  in  two 
regions:  from  0.05  <  x/c  <  0.4  on  the  upper  sur¬ 
face  and  on  the  aft  parts  of  both  surfaces.  In 
the  latter  case  the  shape  of  the  predicted 
pressure  distribution  is,  however,  correct  In  the 
main  and  the  differences  in  level  may  be  ascribed 
to  differences  between  the  measured  and  predicted 
separation  positions.  Using  the  experimentally- 
supported  criterion  (Section  5)  that  a  shape  para¬ 
meter  of  four  identifies  the  point  of  separation 
(ie  the  first  appearance  of  reversed  flow  in  the 
time-mean  velocity  profile),  then  separation  in 
the  experiment  occurs  just  ahead  of  x/c  -  0.80 
whilst  in  the  calculation  it  occurs  Just  ahead  of 
x/c  *  0.90. 

There  are  several  possible  explanations  for  the 
differences  between  theory  and  experiment:  the 
change  in  momentum  thickness  at  the  transition 
trip  could  be  incorrect;  the  tunnel  interference 
induces  a  camber  in  the  flow;  the  three- 
dimensional  nature  of  the  flow  in  the  wind  tunnel, 
and  perhaps  its  slight  unsteadiness,  might  be 
equivalent  to  a  change  in  the  effective  incidence 
of  the  model;  or  there  may  be  significant  defi¬ 
ciencies  in  the  shear  layer  calculations  such  as 
the  omission  of  normal  stress  terms.  The  impor¬ 
tance  of  the  first  three  of  these  possibilities 
has  been  checked  by  appropriate  calculations. 

Increasing  the  momentum  thickness  at  the  trip  by 
0.0003  m  instead  of  0.0002  m  moves  the  separation 
point  forward  to  x/c  ■  0.82  but  the  lift  coef¬ 
ficient  falls  to  1.38  which  is  some  6X  below  the 
experimental  value. 

Using  Fig  1,  the  flow  curvature  and  inclination 
Induced  by  tunnel  Interference  have  been  coun¬ 
teracted  by  adding  an  opposing  camber  to  the  aero¬ 
foil  and  slightly  increasing  its  incidence.  The 
pressure  distribution  about  the  distorted  aerofoil 
at  the  corrected  incidence  has  been  compared  with 
the  uncorrected  pressure  distribution.  Separation 
moved  forward  to  x/c  -  0.86  and  the  consequent 
reductions  in  pressure  near  the  trailing  edge 
contrltated  to  a  slight  (1.4Z)  Increase  In  lift 
coefficient,  but  there  was  no  noticeable  improve¬ 
ment  in  the  pressure  distribution  closer  to  the 
leading  edge. 

A  calculation  at  an  angle  of  incidence  of  13.5* 
(without  an  Induced  camber  correction)  gave  much 
closer  agreement  with  the  experimental  pressure 
distribution  as  shown  In  Fig  19.  The  point  of 
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In  the  calculation  of  the  turbulent  boundary  layer 
the  effect  on  turbulence  structure  of  longitudinal 
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separation  moved  forward  to  x/c  •  0.79  and  the 
lift  coefficient  increased  slightly  to  1.490  which 
is  still  close  to  the  experimental  value.  The 
displacement  thicknesses  from  experiment  and 
calculation  at  angles  of  incidence  of  12. IS*  and 
13.5*  are  compared  in  Fig  20.  The  calculations  at 
the  higher  Incidence  much  more  closely  reflect  the 
trend  of  the  measurements.  This  is  to  be  expected 
in  view  of  the  effect  of  the  incidence  change  on 
the  pressure  distribution.  However  the  measured 
momentum  thicknesses  are  better  represented  by  the 
calculation  at  the  lower  incidence  (Fig  16)  for 
the  rear  part  of  the  wing.  In  this  region  the 
measured  shape  parameter  (H),  which  depends  on 
both  thicknesses,  consequently  falls  between  the 
calculated  results  for  the  two  incidences 
(Fig  21). 

It  is  suggested  that  further  analysis  is  required 
before  the  discrepancy  between  measured  and  calcu¬ 
lated  momentum  thicknesses  can  be  properly 
assessed.  On  the  one  hand  the  effect  of  using  the 
generalised  definition  to  evaluate  the  experimen¬ 
tal  values  needs  to  be  estimated  and  on  the  other 
the  effects  of  Reynolds  normal  stresses  need  to  be 
Included  in  the  calculations. 

8  Concluding  Remarks 

The  flow  near  a  stalling  NACA  4412  aerofoil  has 
been  studied  both  experimentally  and  theoreti¬ 
cally.  The  experimental  results  yield  a  fairly 
detailed  picture  of  the  flow  around  the  aerofoil 
and  in  the  upstream  part  of  its  wake.  The  results 
emphasise  how  the  shear  flow  near  the  rear  of  the 
aerofoil  alters  the  effective  thickness  and  camber 
of  the  disturbance  that  the  aerofoil  Introduces 
into  the  free  stream.  They  also  confirm  that  the 
shear  layers  near  and  beyond  separation  are  thick 
enough  to  require  a  higher  order  theoretical 
treatment  than  traditional  boundary  layer  theory. 

Comparison  of  the  experiment  with  a  semi-inverse, 
viscous-lnvlscld-lnteractlon  theory  shows  that 
lift  coefficients  for  this  aerofoil  can  be  pre¬ 
dicted  very  well  up  to  and  rather  beyond  maximum 
lift.  Predictions  of  detailed  pressure  distribu¬ 
tions  near  maximum  lift  are  slightly  less 
accurate,  though  the  level  of  agreement  with 
experiment  is  very  encouraging. 

Several  possible  reasons  for  the  discrepancies 
between  theory  and  experiment  have  been  advanced 
and,  in  some  cases,  their  likely  Importance  has 
been  tested  by  calculation.  However,  it  is  not 
clear  to  what  extent  the  discrepancies  are  attri¬ 
butable  to  imperfections  In  the  experiment  or  to 
the  need  for  better  representations  of  separating 
shear  layers  in  the  theory.  The  further  develop¬ 
ment  of  this  empirical  aspect  of  the  theory,  if 
Indeed  it  is  significantly  in  error,  will  depend 
on  the  availability  of  suitable  experimental  data, 
such  as  those  presented  here,  for  detailed 
analysis. 
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TABLE  1 

TUNNEL  ROOF  AND  FLOOR  PRESSURE  COEFFICIENTS 


|  *x/c 

1 _ 
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C(VRoof 

( VFloor  j 

-1.0 

-0.079 

+0.106  | 

-0.8 

-0.111 

0.127  j 

|  -0.6 

-0.157 

0.147  | 

|  -0.4 

-0.193 

0.168  j 

-0.2 

-0.226 

0.188  J 
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-0.256 

0.190  | 

1  +0.2 
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0.190  | 

1  0.4 
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0.6 

-0.270 

0.170  | 

0.8 
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0.159  j 

1.0 
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0.131  j 

1.2 

-0.166 

0.107  | 

1.4 

-0.134 

0.089  j 

1.6 

-0.099 

0.040  j 

_ 1 

♦Along  the  tunnel  axis,  here  only. 
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Fig  1  3 


Fig  13  Singularity  model  used  in  panel  method 
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Fig  12  Boundary  conditions  and  corrections  for  EIF 
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Fig  7  First -order  momentum  balance  for 
upper  surface  boundary  layer 
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TABLE  2 

BOUNDARY  LAYER  CHARACTERISTICS 
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- 
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- 
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- 
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•Early  LDA  Profile  Data 


TABLE  3 

HAKE  CHARACTERISTICS 
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Fig  15 
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Fig.  15  Flow  diagram  for  semi  -  inverse  method 


Fig  17  Variation  of  C,  with  <X  for  present  test 
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Development  of  0*/c  on  upper  surface 
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